arXiv:1502.05839v2 [cond-mat.str-el] 5 Nov 2015 


Spiral magnetism in the single-band Hubbard model: the Hartree—Fock and 

slave-boson approaches 

P.A. IgoshcQ 

Institute of Metal Physics, 620990, Kovalevskaya str. 18, Ekaterinburg , Russia and 
Ural Federal University, 620002 Ekaterinburg, Russia 

M.A. Timirgazin, V.F. Gilmutdinov, and A.K. Arzhnikov 
Physical-Technical Institute, 426000, Kirov str. 132, Izhevsk, Russia 

V.Yu. Irkhin 

Institute of Metal Physics, 620990, Kovalevskaya str. 18, Ekaterinburg, Russia 

(Dated: November 6, 2015) 

The ground-state magnetic phase diagram is investigated within the single-band Hubbard model 
for square and different cubic lattices. The results of employing the generalized non-correlated mean- 
held (Hartree-Fock) approximation and generalized slave-boson approach by Kotliar and Rucken- 
stein with correlation effects included are compared. We take into account commensurate ferromag¬ 
netic, antiferromagnetic, and incommensurate (spiral) magnetic phases, as well as phase separation 
into magnetic phases of different types, which was often lacking in previous investigations. It is 
found that the spiral states and especially ferromagnetism are generally strongly suppressed up to 
non-realistically large Hubbard U by the correlation effects if nesting is absent and van Hove sin¬ 
gularities are well away from the paramagnetic phase Fermi level. The magnetic phase separation 
plays an important role in the formation of magnetic states, the corresponding phase regions being 
especially wide in the vicinity of half-filling. The details of non-collinear and collinear magnetic 
ordering for different cubic lattices are discussed. 

PACS numbers: 75.10.Lp, 71.10.Fd, 71.10.Hf, 75.30.Fv 


I. INTRODUCTION 

Magnetic properties of strongly correlated transition-metal compounds and their relation to doping, lattice geometry, 
band structure and interaction parameters are still being extensively investigated^. In particular, the details of 
magnetic order in the ground state remain to be examined both theoretically and experimentally. During recent 
decades, the two-dimensional (2D) case closely related to the problem of high-temperature superconductivity in 
cuprates and iron arsenides has been intensively investigated theoretically. However, the more widespread case of 
three-dimensional compounds is not yet sufficiently understood. In general, the following factors should be taken 
into account to state an appropriate physical picture of the magnetism formation: (i) spiral magnetism, (ii) phase 
separation, (iii) electron correlation impact, (iv) multiorbital effects. 

The single-band Hubbard model provides a sufficient basis for an analysis of complex magnetic properties. In 
the case of a large Coulomb interaction parameter the ground state of single-band model in the half-filled case for 
bipartite lattices is a Neel antiferromagnetic (AFM) insulator. The types of the antiferromagnetic state instability in 
the presence of doping or the finite integral of the electron transfer between the second neighbors are still incompletely 
revealed. According to the classical work by NagaokaA, on addition of one excess charge carrier the ground state on 
the bipartite lattice becomes saturated ferromagnetic (FM). This statement can be also considered as a reasonable 
hypothesis in the case of small finite doping^ - —. At the same time, for the non-bipartite (e.g., fee and hep) lattices the 
ground state can be more complicated^. For finite doping, the ferromagnetic state can become non-saturated^—. 

Scenarios of the possible doping-induced magnetic ordering include the phase separation (PS) of different types: 
into the ferromagnetic and antiferromagnetic phases^ or the phase of superconducting electron liquid and the Neel 
antiferromagnetic phased. An alternative scenario is the formation of a spiral (incommensurate) magnetic state. 
It was considered within different approaches: the analysis of the momentum dependence of the generalized static 
magnetic susceptibility for the bare spectrunJA, the Hartree-Fock approximation (HFA) (with small and moderate 
U/W values considered, where U is the Coulomb repulsion parameter and W is the bandwidth) 12 : 13 and the t — J 
model (large U/W values)^ -15 . It should be stressed that within simple HFA-like approaches the paramagnetic (PM) 
state is, as a rule, less favorable than the magnetic ones, which seems unreasonable because the PM phase is dominant 
in the ground state for the majority of transition metal compounds. Another issue is investigating competition of the 
collinear, incommensurate and commensurate (FM and AFM) phases. 

Non-uniform phase-separated states are observed in both 2D and 3D magnetic system s 16 : 17 . The experimental 
observation of antiferromagnetic or spiral magnetic structures makes the classical problem of theoretical description 
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of the formation of different types of magnetic order very urgent. Spiral structures are observed in both two- and 
three-dimensional compounds: in iron-based high-temperature superconductors!! and fee iron precipitates!!, MnSi 
systems!!, itinerant compounds Y 2 Fei 7 and Lu 2 Fei 7 22 ' 23 where helical magnetism is induced by pressure. Incommen¬ 
surate structures are found in AFM chromium!!. In doped cuprates such states are observed as the dynamic magnetic 
order—. Besides, considerably enhanced incommensurate magnetic fluctuations are observed in strontium ruthenates 
at low temperatures!! (see discussion in Refsd 3 ' 26 ' 27 ). 

Although spiral states were widely theoretically studied, the relation between the type of magnetic order and the 
model parameters has not been still clearly determined. In addition, there is a question about the possibility of 
the inhomogeneous state formation in a broad sense: either via magnetic phase separation or via so-called stripes. 
The latter structures, representing a mixed charge-spin order, are widely discussed for 2D high-T c superconducting 
systems!! - —. 

A detailed study of the magnetic phase diagram of the 2D Hubbard model taking into account the electron transfer 
integral only between the nearest neighbors, t , within HFA shows that the spiral magnetic states occur in a wide range 
of parameters, especially at moderate values of U < W 31 ' 32 . It was shown in Ref. [H that inclusion of the electron 
transfer between the second neighbors, t' ^ 0, in the Hamiltonian changes considerably the ground-state magnetic 
phase diagram. It was established that the inhomogeneous stripe structures can be more energetically favorable than 
the spiral states at t' ^ 0 in the slave-boson approximation!!, as well as at t' = 0 and U < W in the Hartree-Fock 
approximation— However, this conclusion is somewhat devaluated by disregarding the intersite Coulomb interaction 
which can considerably increase the energy of the stripes. 

Last decades, the dynamical mean-field theory (DMFT) has been intensively developed and became a powerful tool 
for the solution of quantum many-body problems— Being well justified in the limit of large space dimensionality d, 
it takes into account local correlations exactly, which can, in principle, provide some unbiased solution of the problem 
of ground magnetic state in the Hubbard model. In Refs. 13611371 a DMFT generalization to spiral magnetic states was 
proposed. A typical HFA pattern of phase transitions between commensurate and incommensurate magnetic states 
with doping was reproduced for the square lattice. A novel so-called inhomogeneous DMFT based approach was 
proposed in Ref. [38) to study magnetic state formation for a rather large cluster. The collinear stripe-like spin-density 
wave states were found at finite doping and the magnetic phase diagram for the square lattice was constructed. 

Another convenient approach to study the magnetic order formation with account of correlations is the slave-boson 
approximation (SBA) by Kotliar and Ruckenstein!! within the Hubbard model. It yields a simple (in computational 
sense) method taking into account the crucial difference between the single and double site states on the mean-field level 
and allows a thermodynamical description of the magnetic phase transition physics. In the saddle point approximation, 
this method is qualitatively similar to the Gutzwiller approximation's. The ground state energy obtained within the 
slave-boson method is in a good agreement with the quantum Monte-Carlo and exact diagonalization calculations— 
This justifies the use of the slave-boson method in the cases where there are no strong fluctuations of the order 
parameter in the system. In the opposite case of the continuous (second-order) phase transition, where strong 
fluctuations occur, the nontrivial many-particle renormalizations of the Green’s functions can be importan t 27 ' 42 . 

The effect of electron correlations on the stability of spiral magnetic states was considered in— using the SBA by 
Kotliar and Ruckenstein. The phase diagram of the Hubbard model was built in the nearest-neighbor approximation 
( t' = 0). Derivation of the static magnetic susceptibility in SBA— allowed one to obtain a criterion for the PM state 
instability with respect to the second-order transition into an incommensurate magnetic state. This generalizes the 
analogous criterion obtained within the random phase approximation 44 to the strongly correlated states. Using the 
generalized criterion, the finite t' case was studied and a considerable tendency to FM ordering at hole doping and 
large t'/t values was found!!. However, the phase transitions between magnetically ordered states, as well as all 
first-order transitions cannot be studied within this approach. 

The slave-boson approach can be generalized to take into account the multi-orbital local interaction part of the 
Hamiltonian with an indirect connection between the corresponding eigenstates and slave-bosons states!!. The method 
proposed (the so-called spin-rotation-invariant SBA) yields the effective momentum dependence of the Green’s func¬ 
tion residue which is missing in the traditional slave-boson approaches. Recently, a successful combination of the 
slave-boson approximation and the local density approximation (LDA) was used for the first-principle investigations 
of real compound s 47 ' 48 as a numerically inexpensive way for qualitatively correct treatment of the electron correlation 
in the electron spectrum beyond LDA. Further investigations generalized the method for the study of superconductiv¬ 
ity!!, the physical picture being similar to that obtained within multiorbital LDA+DMFT approach. The excitation 
spectrum in the 3D case differs essentially from the 2D case, and finite-temperature behavior can be hardly described 
by SB method. However, in the ground state the SB approach seems to be useful. To generalize the SBA to the 
case of finite temperatures the Gaussian fluctuations can be included, e.g. in the framework of statistically consistent 
Gutzwiller approximation (SGA)— 

In the present paper we apply both the non-correlated mean-field (Hartree-Fock) approximation and the slave- 
boson approximation to the magnetic ordering of different types (ferromagnetic, antiferromagnetic and spiral) within 
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2D and 3D Hubbard model and compare the results. In Section mi we present the model and approximations used. 
The main physical factors which favor the spiral structures formation are discussed in Section 1IIII Sections IIVI and [V] 
are devoted to the discussion of the results obtained for the 2D and 3D case, respectively. Conclusions are presented 
in Section |VT] 


II. MODEL AND SLAVE-BOSON APPROACH 

We consider the Hubbard model 


^ ^ ' tijSfja-i c ia Cj a , ~b U ^ (1) 

ijcrcr' i 

where the matrix elements of the electron transfer are tij = —t for the nearest neighbors and t' for the next-nearest 
neighbors (we assume t > 0), cj a , c ia are the electron creation and annihilation operators,respectively, = c\ a c ia , i 
is the site number, a is the spin projection. 

The local spin space rotation around x axis, matching different site magnetization vectors along, say, z axis, by the 
angle QR^ (where Q is a spiral wave vector, Ri is the site position) is applied for the consideration of magnetic spirals. 
This maps the spiral magnetic state into an effective ferromagnetic one, but the hopping term in the Hamiltonian 
becomes non-diagonal with respect to index a: tij5 ac ,' —> tf? in Eq. ©■ 

^ +U J2 n 4 n 4> ( 2 ) 

ijcrcr' i 


where tf? = exp[iQ(R, — R j)^ x ]<ra'tij- 

Further Hartree-Fock treatment of the many-particle Coulomb interaction term replaces it by the one-electron 
interaction term with some effective field U{riis) which however mixes the averaged contributions from singly and 
doubly occupied states (( Ui a ) is the statistical averaged electron density at site i and spin projection a). This is not 
correct even qualitatively, especially at large U. 

Since doubly occupied site states are energetically unfavorable, a simple way of taking into account the correlation 
effects is to introduce the slave-boson operators ei(e |) , pi a (pl a ), di{d\) for empty, singly and doubly occupied states, 
respectively^. The transitions between the site states originating from intersite electron transfer can be formulated 
in alternative terms, but simultaneously with the conventional Slater determinant based formalism (related to the 
electron creation/annihilation operators). Conceptually this is close to the Hubbard’s X-operator formalism^i, 
where, however, the site transition X-operators are introduced instead of the original one-electron operators. We use 
the action formalism and pass to the extended Hamiltonian doing in the fermion-boson Fock space, the Coulomb 
interaction term becoming diagonal with respect to the boson operators (see method details in Ref. [39): 

Hext = J2 C L C j*' Z L Z j<r' + U J2 d * d i’ ( 3 ) 

ijcrcr' i 


where 


z ia = (1 - d\di -pI^) 1/2 {e\p ia +p\ s d i ){ 1 - eje* -pt_ Pi _) V 2 . (4) 

Equation ((3]) is equivalent to Eq. m provided that the unphysical site states (which have no counterparts for the 
original electron system) are excluded (Eq. d5j) below) and the electronic and bosonic subsystems have a strict mutual 
correspondence (Eq. © below). The physical subsector of the domain of operator H ex t is restricted by the following 
constraints: 


e l e i + Yl P ^Pia + 4 d i = 1 > ( 5 ) 


4 d i+Pi*Piv = c L c iv (6) 

The presence of the constraints can be taken into account within the functional integral formalism through the 
Lagrange multipliers (the local bosonic “chemical potential” ry for Eq. © and the on-site electron-boson “force” A i G 
for Eq. ©) introduced into the action. To this point the transformation considered has been rigorous, but further exact 
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functional treatment of both bosonic and fermionic fields is hardly possible. A reasonable physical picture (different 
from the Hartree-Fock one) at arbitrary U/t can be obtained within the saddle-point approximation for the bosonic 
fields and Lagrange multipliers: for the action S generated by Eq. © we replace the bosonic fields by their extreme real 
values which are assumed to be site- and time-independent: e \, e, : —> e; p\ a , p ia . —> p a ; d\, d 1 —t d; pi — > p, Xt a —> X a . 
Explicit equations for the extreme values are presented below. Note that within this approximation the numbers 
e,p a ,d have a simple meaning: these are the probability amplitudes of any site’s being in the corresponding many- 
electron state. 

The fermionic part of the action can be produced by the Hamiltonian 

Tdf = 'y ( {ZqZg’t^ + SijS aa ' X a )c\ a .Cj a i , ( 7 ) 

cra'ij 

the corresponding eigenvalues (antiferromagnetic subbands) yielding the renormalization of the electron spectrum: 

-E«=±i(k) = (1/2) ((z 2 + zf)e £ + Af + Aj,) + s^/T>k, ( 8 ) 

where 

-Dk = (1/4) ((z 2 — zf)e£ + A-|- — A 4 .) + (e k ztA |.) 2 (9) 

and 

e k = (!/ 2 )(^k+q/2 ± ^k-q/2), (10) 


= exp(ik(R ?; - 


( 11 ) 


is the bare electron spectrum. The subband narrowing is the c-number function of extremal bosonic fields 

4 = (!- d 2 -pl){e P<J +p s d) 2 ( 1 - e 2 -pi). 


( 12 ) 


The impact of the averaged site states on electron states manifests itself in two types of bare spectrum renormaliza¬ 
tion: (i) narrowing of the bare spectrum (similar to the Hubbard I approximation^) which is specified by the factor 
z 2 < 1, (ii) the additional energy shift X a . Note that within HFA 

A ct A« f = Un/2 - Uma/2, (13) 

where 

n = J2( c L c ia)> m = Yl a ( c L c ia) ( 14 ) 


are the electron density and the site magnetization. Both z 2 and A CT are essentially spin dependent, which allows 
the magnetic states formation to be studied. However, in the slave-boson formalism, in contrast with HFA, the one- 
electron energy shift A CT cannot be expressed in terms of n and m only and should be specified separately, see Eq. 
(fT71) below. 

Then we can write the saddle-point (mean-held) equations: 


n 


m 

Xa 


ks 

1 x—' et(zi 

-y>/(£,( k»-SU: 

ks 


"f) + At ~ Aj. 

VDk 



_d/Pa_\ 

P 2 a+d 2 ) 


/{epa) , 


(15) 

(16) 
(17) 


where we have introduced v = ed — p^p\,^ f(E) = (exp((E — p)/T ) + 1) 1 is the Fermi function, p is the chemical 
potential, N is the lattice site number. We have introduced the quantity 


= 


epa + Pad 1 

(e 2 +p 2 a){p 2 a +d 2 ) N 




ks 


dE s ( k) 
d(z 2 ) ' 


( 18 ) 
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which is an analogue of the Harris Lange shift^, and the derivative dE s (k)/d(z 2 ) corresponds to the formal derivative 
of Eq. © with respect to z 2 . 

The saddle-point values of bosonic variables should be obtained as a solution of the self-consistent equations: 


'£ 


2 d 2 +p 2 +p 2 = 

n, 

(19) 

2 2 
Pt~Pi = 

m, 

(20) 

e 2 +p 2 +p 2 + d 2 = 

1, 

(21) 

(— + = 
\epa Pad J 

U. 

(22) 


We consider now formal relation between Hartree-Fock and slave-boson approximations. Note that the equations 
for HFA can be represented in a form similar to SBA, since the former one treats the double occupancy as 

d 2 = (d\di) = (n^ni) ( n t )(ni ), (23) 

which can be replaced in virtue of Eqs. m an d h2qd, 

d 2 = ( P 2 + d 2 )(p 2 + d 2 ), (24) 

which implies 


d 2 (e 2 ) 




(^) 


2 


2 2 


(25) 


where the signs are chosen to satisfy d 2 < e 2 for n < 1 and vice versa. This replacement of Eq. (E21) by Eq. (E3 
allows to formulate HFA approximation in terms of bosonic quantities. Now one can see that the HFA differs from 
SBA because the former neglects the different nature of the singly and doubly occupied states forming the electron 
spectrum. The rough ansatz (1231) does not include the bandwidth narrowing ( z 2 factor) caused by Pauli principle, 
which is taken into account within SBA: Eq. (l?5l) implies that v = ed — p^p\_ = 0 and, from Eq. (fl2l) . z2 = z? = 1. 
However, Eq. (l22l) means that v = 0 implies X a = 0 ^ A^ F , which manifests the inconsistency of Eqs. m and (l23l) 
within HFA. A concrete example for the concentration dependences of bosonic fields defined above within HFA and 
SBA is presented in Sect. IV PI 

The saddle-point expression for the spiral state thermodynamic potential (per site) f2(/i, Q) has the form 


H( m , Q) = Ud 2 -J2 A a(p\ + d 2 ) + n f (p, zl, A ct ), (26) 

a 

where ttf(/j,,z 2 , X a ) is the standard potential for the non-interacting electron system in the field Act with a narrowed 
(by the factor z 2 ) spectrum 


(4, Act) = ~(T/N) J2 log (1 + exp(—(-EA(k) - p)/T )). (27) 

k s 

The resulting wave vector is determined by the minimization of Q over various spiral states (the spiral wave vector 
Q) at fixed p in the limit T —> 0. The minimization of fl(Q) with respect to Q was performed numerically with Q 
running in high-symmetry directions of the Brillouin zone (we choose the lattice spacing as unity): the (Q,Q), ( Q , n) 
and (0, Q) directions were considered for the square lattice; the ( Q , Q, Q), ( Q , Q, i r), (Q, n, n), (0, 0, Q), (0, Q , Q) and 
(0, <5,7t) directions were taken into account for cubic lattices (for bcc and fee lattices 7r is replaced by 27r due to 
halving the lattice spacing). The step of changing Q was 0.027T for the square lattice and 0.057T for cubic lattices. The 
number of k-points in the Brillouin zone during integration was 400 for the square lattice and 70 for cubic lattices. 
The accuracy settings were adjusted, if necessary. Since f 1 in the ground state actually depends on the chemical 
potential p, as a parameter, we can determine the dependence of the magnetic structure on p by a procedure taking 
into account the possibility of the phase separation-^, avoiding the tedious application of Maxwell rule in the negative 
compressibility regions. 

To illustrate the non-uniform state formation we present an example of the density dependence of the chemical 
potential, calculated by commonly used method and the present method (see also Ref. 113! ) in the vicinity of half-filling. 
One can see that the former method yields a thermodynamical unstable dependence p{n) which should be corrected 
by additional applying the Maxwell rule. 
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FIG. 1: (Color online) The dependence of the chemical potential /j,, magnetization m and Q x vs electonic density for simple 
cubic lattice, calculated within SBA, t' = 0, U = 14t, Q = (Qx, tt, n). Solid line presents the result with account of PS, dashed 
line is the standard result of Ref. |u]- The choice of the wave vector form provides the minimum of the total energy (see Fig. [4] 
below). 


III. FERMI SURFACE NESTING AND VAN HOVE SINGULARITIES 


Before proceeding to the results obtained we need to mention the main physical factors that are important for the 
incommensurate magnetic states formation starting from paramagnetic phase. There are two well known mechanisms 
leading to stabilization of incommensurate states: (i) Fermi surface nesting and (ii) van Hove scenario. 

(i) Fermi surface nesting is the presence of large segments of a Fermi surface that can be connected to other segments 
via a nesting vector Q in such a way that occupied states move to the empty states region and vice versa. When the 
nesting condition et = — £k+Q is realized for many k-points the magnetic susceptibility function 


k 


/( e k) - /(ek+Q) 
£k+Q — Ck 


(28) 


is enhanced or has a singularity leading to stabilization of the incommensurate magnetic state with wave vector Q^. 
One can mention that the nesting condition means = 0 in our formalism (see (1101) ). This results in cancellation 
of the symmetric terms in spectra 2?(k) (|8|). and hence the ^-dependence of E( k) becomes appreciably suppressed. 
Thus we come to an important conclusion that the Fermi surface nesting makes the system to some extent insensitive 
to the correlation effects. In the FM case Q = 0 and the nesting condition reduces to ek = 0. This means that the 
larger is the density of states at a given /x, the weaker is the influence of electron correlations on the FM state. 

(ii) Another mechanism leading to the incommensurate magnetic structure formation was proposed in Ref. I&H . 
The authors have shown that the presence of saddle points in the electron spectrum results in the divergence of 
susceptibility (051) if Q is equal to a vector connecting the saddle points to each other. The result was obtained for 
the 2D case but one can expect that it holds true for the 3D case as well with a possible change of the divergence 
into a weak singularity. This mechanism demonstrates the importance of van Hove singularities (given by the saddle 
points) for magnetism. 

Thus, when analyzing the magnetic phase diagrams, special attention should be paid to the density of states, nesting 
properties of the Fermi surface, and van Hove singularities. 


IV. RESULTS: SQUARE LATTICE 

A. t' = 0 

We calculated the phase diagrams of the ground state at different t'/t values within the Hartree-Fock approximation 
in Ref. [HI; the results for t' = 0 are presented in Fig. (5Ji. An analogous phase diagram obtained within SBA is shown 
in Fig. (we restrict ourselves to the case 0 < n < 1 due to the particle-hole symmetry, see justification below). 
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Since SBA takes into account the impact of many-electron states onto the one-electron spectrum on the static level, 
we can treat these SBA corrections as correlation effects beyond HFA. Note that for both these approaches the phase 
transitions between different magnetic states are generally first-order transitions, which leads to considerable regions 
of phase separation. This also applies to the transition from AFM to spiral state at arbitrary small doping, which 
is of the first order due to that the insulator finite-gap state goes abruptly to the metallic phase (see also22). The 
account of PS substantially distinguishes the phase diagram obtained within SBA from that presented in Ref. Ell 
where PS was disregarded. The separation regions between the AFM phase and the spiral magnetic states (parallel, 
Q = (Q, 7 r), and diagonal, Q = (Q, Q)) are rather wide; the corresponding regions of pure spiral states are narrowed. 
In particular, this refers to the diagonal phase, the existence of which becomes possible only at U > lit. The phase 
transition between the PM and spiral magnetic states is of the second order. The magnitude of number Q in the 
regions of spiral states is variable but qualitatively it behaves as shown in Fig. [T) the larger is the doping, the smaller 
is the value of Q. This behavior of Q is typical for all the lattices considered. 


Hartree-Fock Slave boson 



FIG. 2: (Color online) Ground state magnetic phase diagram of the Hubbard model for the square lattice with t' = 0 at 
n < 1 within (a) HFA^ and (b) SBA. The spiral phases are denoted according to the form of their wave vector. Filling 
shows the phase separation regions. Bold (blue) lines denote the second-order phase transitions. Solid (red) lines correspond 
to the boundaries between the regions of the homogeneous phase and phase separation. Dashed (red) lines separate the PS 
regions corresponding to different phase pairs. The second-order phase transition boundary produced by the Stoner criterion 
U c = 1/p(Ef) (where p(Ep) is the density of states at the Fermi level in the paramagnetic phase), separating the ferromagnetic 
and paramagnetic regions, is shown by dash-dotted (violet) line. Qafm = (n,ir). 

One can see that the electron correlations lead to a noticeable suppression of the magnetically ordered states in 
comparison with HFA: the corresponding concentration intervals in the phase diagram decrease strongly, and the 
variety of the spiral states disappears. Besides, in the slave-boson approximation there occurs a wide region of PM 
state. The FM state covering a considerable part of the diagram within HFA is shifted to the region of large values 
U > 60t. This is a manifestation of the fact that SBA yields a convenient scheme of interpolation between the 
small U limit (coinciding with HFA) and the large U limit (d = 0 as the limit consequence of Eq. (1221) 1 in the case 
n < 1; whereas HFA yields only saturated ferromagnetism at large U for finite doping. The suppression of FM state 
reproduces the result obtained in Ref. El) and is in good agreement with the variational study of the stability of the 
saturated FM phased. The region of separation of the AFM and spiral phases is narrowed by about one half (with 
respect to the variable n). 

According to our calculations, even an unrestricted increase of U does not make the magnetically ordered states 
stable far from half-filling: at U —> oo, there are no spiral magnetic solutions of the slave-boson method equations 
at 11 — rz.| > 0.63. At the same time, the saturated FM solution becomes more favorable than the spiral ones only 
at |1 — n| < 0.3. This improves the result proposed in Ref. where spiral phases were not considered; similar 

phase diagram was obtained in Ref. ED but PS was disregarded. Thus, at large U values far from half-filling the 
spiral magnetic state replaces the saturated FM one. In contrast to Ref. 0, in our approach there exist unsaturated 
FM solutions, but they are always energetically unfavorable in comparison with the saturated FM or spiral magnetic 
states. 

In order to compare our results with the Stoner criterion and analyze the van Hove singularities, the inverse density 
of states (DOS) at the Fermi level is depicted in Fig. [2] and all diagrams below. There is a saddle point in the electron 
spectrum of the square lattice k = ( 0 , 7 r) which produces a divergent van Hove singularity of DOS. For t' = 0 this 













singularity is at the band center (we introduce the notation for the density corresponding to a van Hove singularity 
in the paramagnetic spectrum n v hi = 1) and coincides with the case of perfect nesting. Both these peculiarities 
determine the weakness of correlation effects in the vicinity of half-filling, cf. Fig. [l] which is reflected in qualitative 
accordance of the HFA and SB A results close to ?r V Hi- 


B. t' = 0.2 1 

Up to the end of this Section we present our results for the square lattice for the nonzero value of t' at which the 
particle-hole symmetry breaks down: of great interest is the difference between the properties of the hole-doped (n < 1) 
and electron-doped (n > 1) systems^. The consideration of the negative ratio t'/t can be reduced, as a consequence of 
lattice bipartiteness, to that of the positive one by the particle-hole transformation c ia —f (—1 ) 1 c| cr , t' — > —t', n —> 2—n, 
where the site numbering provides different parity for neighboring sites. We present the phase diagram for t’ = 0.2 t 
within HFA (Fig. [3ji) and SB A (Fig. [3 }d). In comparison with the case of t' = 0, for n < 1 the diagonal phase shifts to 
the region of much smaller U/t values and the parallel phase becomes more extended over the concentration parameter. 
The physical cause is the shift of the density value n V H (for t' = 0.2 t n V Hi ~ 0.83) from the half-filling (n = 1): at 
t' ^ 0, tt-vH and the half-filling determine two different trends to the corresponding magnetic orderings. We see that 
in the hole-doped half of the phase diagram the correlation effects lead only to the quantitative renormalization of 
the phase boundaries. 


Hartree-Fock Slave boson 



FIG. 3: (Color online) Ground state magnetic phase diagram of the Hubbard model for the square lattice with t' = 0.2t within 
(a) HFA^; (b) SBA. Bold dotted (red) lines denote the first order phase transitions if the phase separation region is narrow. 
The notations are analogous to those in Fig. [2] 

In the case of n > 1, the correlation effects are more significant: all homogeneous spiral states disappear, except for 
a narrow region of the parallel phase. Such a suppression of magnetism can be explained by the fact that in this case 
all peculiarities (imperfect nesting) come from the point n = 1. Far from half-filling, no spiral magnetism is possible 
and the saturated ferromagnetism is energetically unfavorable in comparison with the PM phase at any large U value, 
when correlations are taken into account. Thus, the particle-hole asymmetry is enhanced considerably as compared 
with HFA. 

These results can be used for qualitatively explaining the magnetic properties of layered high-temperature super¬ 
conducting perovskites, for which the fitting of the angular resolved photoemission spectroscopy spectra to those of 
the bare model yields t' ~ 0.2 t, se o 13 : 54 and references therein. The results are in agreement with the experimental 
data on the magnetic structure of the hole-doped compound La 2 - p Sr p Cu 04 which has a close value of the t'/t pa¬ 
rameter—At the same time, for the high-temperature superconducting compound Nd 2 - rc Ce a; Cu 04 , in which the 
charge carriers are electrons, the homogeneous commensurate AFM ordering is stable up to x = 0.14^ in agreement 
with our results for n > 1 (see Fig. 0. 
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V. RESULTS: 3D LATTICES 

Now we present the results for three-dimensional lattices: simple cubic (sc), body-centered cubic (bcc), and face- 
centered cubic (fee) lattices. We restrict ourselves to the collinear and spiral structures with magnetic cell being not 
increased, so that.,e.g., type-III antiferromagnetic collinear structure with a doubled cell^, is beyond our consideration; 
such structures will be treated elsewhere. 


A. Simple cubic lattice 

1. t' = 0 


Hartree-Fock Slave boson 



FIG. 4: (Color online) The same as in Fig. [2] for the simple cubic lattice Hubbard model with t' = 0. Qafm = (7r,7r,7r). 

The physical picture is very similar to that for the square lattice (see Fig. 0 for a comparison of the results of 
HFA (a) and SB A (b) approaches): the density value n = 1 corresponding to the perfect AFM nesting peculiarity 
(but not to a van Hove singularity) retains its crucial role. The spiral magnetic phase with Q = (Q,7t) of the square 
lattice is replaced by a spiral with Q = (Q,7r,7r). The electron spectrum of the sc lattice has two saddle k-points: 
the X and M points of the Brillouin zone which give two van Hove singularities of DOS at the electron concentrations 
n V Hi ~ 0.425 and n V H2 ~ 1.575, respectively. The vector connecting two X-points with coordinates (0,0,7r) and 
(0, —7T, 0) determines the tendency to Q = (0, n, n) state by the van Hove scenario described in Section Hill In contrast 
to the 2D case, the van Hove singularities are bounded and do not much affect the SBA diagram — the correlation 
effects strongly suppress magnetism at n being far from half-filling, including the regions close to n V H’s. 


2. t' = 0.3 1 

In Figs. and 0> are presented the respective HFA and SBA results for the case of positive non-zero t'/t, when 
these tendencies come out in different parameter regions on the phase diagram (we consider the case of t’ = 0.3f, the 
consideration of the negative ratio t'/t is omitted for the same reason as in Section HVl) . 

This case is characterized by four van Hove singularities (all being bounded): n V Hi ~ 0.23 (k = (0, a, a), cos a = 
t/2t' — 1); n V H2 ~ 0.26 (k = (/?, /3, /3), cos/3 = t/4t'); n V H3 ~ 0.368 (T point); n V H4 ~ 1-41 (M point). The saddle points 
which give the ’vHl’ and ’vH2’ singularities do not correspond to special high-symmetry points of the Brillouin zone 
and hence the van Hove scenario described in Sec. IIIII favors the formation of incommensurate vertical (Q = (0,0, Q)) 
and diagonal (Q = (Q,Q, Q)) states. The reduced symmetry of the saddle points leads to the removal of their 
degeneracy and, as a result, to an increase of the DOS. The large DOS value stabilizes the FM phase for a rather 
small U/t~ 1.5 and protects it from suppression by the electron correlations (see the SBA diagram). Also the spiral 
states turn out to be stable to correlation effects, which can be explained by the Fermi surface nesting formation. 
The Fermi surface section for the case where the paramagnetic state density n = n V Hi is depicted in Fig. 0 It 
consists of six cylindrical figures with regular squares at their bases. These squares are connected by the nesting 
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vector q = (0, 0,q),q = 2 arccos(t/2t' — 1) « 0.2677T. In accordance with the nesting concept this leads to magnetic 
instability with respect to the formation of incommensurate structure with q = Q. A comparison with the calculated 
phase diagram (Fig. 0 shows an agreement in the Q value with the (0, 0, Q) phase found near n V Hi- 


k z 

-n n 



FIG. 5: (Color online) Fermi surface section in a half of Brillouin zone for sc lattice, t' = 0.3f, n = n V Hi- q is the nesting 
vector, see text. 

As well as in the 2D case, we found a strong hole-particle asymmetry: for n < la variety of magnetic phases appear 
at not large U/t. For instance, at a realistic value U = At we found the AFM phase in the vicinity of half-filling. As 
the density lowers, it is replaced by a phase with Q = (Q, and a small fracture of (Q,Q,tt) phase; on further 
lowering of the density a wide non-magnetic region is observed, and at low density the system enters into the region 
of competition between the FM, ( Q , Q, Q) and (0, 0, Q) phases. 


Hartree-Fock Slave boson 



FIG. 6: (Color online) The same as in Fig. [3] for the simple cubic lattice Hubbard model with t' = 0.3 1 (a) HFA results; (b) 
SBA results. Qafm = 

To illustrate the properties of phase separation into Neel antiferromagnetic at n = 1 and spiral (ferromagnetic) 
phase (n ^ 1) we present their U dependence (z%, m, Q and n ) along the boundaries of PS region (see Fig. 0 for 
both these phases participating in PS within SBA in Fig. 0 As U increases both m and (z|) for n < 1 [n > 1) 
tends to 1, on the other hand, z? ( z |) for n < 1 (n > 1) is strongly reduced. The density n at the boundary (n ^ 1) 
appears to be only weakly [/-dependent. These tendencies are much stronger for the AFM phase at n = 1, which 
implies that correlation effects are much more pronounced at n / 1 than in half-filled case n = 1 where z a are close 
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to unity. The magnetization jump between these two phases turns out to be almost [/- independent, but considerably 
different for n < 1 (|Am| ~ 0.2) and n > 1 (|Am| ~ 0.4). 




FIG. 7: (Color online) Ground state characteristics, z%, m, Q x /tt and n, of the phases participating in PS in the vicinity of 
half-filling for the sc lattice at t' = 0.3t as a function of U: (a) n < 1 and (b) n > 1. Solid lines in (a) and (b) correspond to 
left and right boundary line of PS region, respectively; dashed lines present ttzaf = m and z \ F = z% = zf for the AFM case 
at n = 1. Line breaks in (a) originate from the first order transitions, being connected with the jump of the spiral state wave 
vector Q. The magnetization jump ttiaf — m is given by the difference of dashed and solid lines. 


B. Body-centered cubic lattice 

1. t' = 0 

This lattice is characterized by the density of states with a strong log 2 -singularity at the band center for t' = 0, which 
corresponds to a whole line of van Hove singularities in the Brillouin zone (k = (Q,7 t,7t), 0 < Q < it , D direction). 
In Fig. [8] we present the results for the bcc lattice Hubbard model for both HFA (a) and SBA (b) at 0 < n < 1 
(the particle-hole symmetry holds). The results are somewhat similar to those for the square lattice (logarithmic 
singularity). One can see that the correlation effects strongly but only quantitatively renormalize the phase boundary 
lines in a rather wide vicinity of half-filling, whereas well away from half-filling they fully destroy magnetic ordering. 
The van Hove singularity corresponding to n = n v h = 1, results in this case in especially wide regions of PS and 
magnetic phases. This is a manifestation of the particularly strong bcc lattice van Hove singularity that merges with 
the half-filling peculiarity related to the perfect nesting. The diagonal magnetic order (Q = ( Q , Q, Q )) dominates, 
and at moderate Ujt this is the only existing order. In general the magnetic phase distribution is similar to the case 
of square lattice (see Fig. E). 


2. t' = 0.3f 

We also present the results for the bcc lattice with t' = 0.3 t obtained within HFA (Fig. [9b) and SBA (Fig. 
[2b)- In this case we have 3 bounded van Hove singularities: n v ni ~ 0.45 (P point), n V H 2 ~ 0.494 (k = (a, a, a), 
cos a/2 = t'/t), n V H 3 ~ 1-04 (N point). The variety of magnetic phases which exist far away from half-filling in HFA, 
vanishes. The results are considerably different from those for the square and sc lattices. Despite the presence of van 
Hove singularities at n < 1 the FM ordering is strongly suppressed in favor of spiral ordering. At moderate U > 12i 
the ground state FM ordering appears to be non-saturated (1 > m/n ^ 0.5), and the saturated ferromagnetism occurs 
only at U > 24 1. For hole doping at rather large U/t > 10 the FM phase and the (satellite) Q = (Q, Q,Q), (0,0, Q) 
phases (they both can be considered as modulated FM phases) survive within SBA. At smaller U/t the diagonal 
(Q, Q , Q ) phase becomes favorable in a wide region. 

The AFM phase region is wide (especially in SBA) even away from half-filling for n > 1. The correlation effects 
stabilize pure AFM ordering with wave vector Q = (0, 0, 27 t) up to large U/t in the vicinity of half-filling (1 < n < 1.4). 
This is explained by the nesting on the one hand and by the ’vH3’ singularity on the other hand. Q = (0,0,27r) 
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Hartree-Fock Slave boson 



FIG. 8: (Color online) The same as in Fig. [2] for the bcc lattice Hubbard model with t' = 0, Qafm = (0, 0, 2-7t). 

Hartree-Fock Slave boson 




FIG. 9: (Color online) The same as in Fig. [3] for the bcc lattice Hubbard model with i! = 0.3t. (a) HFA results; (b) SBA 
results. Qafm = (0, 0, 27 r). 


corresponds to the vector connecting (0,7r,7r) and (0,7r, — n) N points of the Brillouin zone. For n > 1.4 the AFM 
state changes to the PM phase through a second order transition. The magnetic phase Q = ( 7 r, 7r, 7 r) that occurs 
within HFA, as well as the corresponding PS region, turn out to be fully destroyed by the correlation effects. Lowest 
critical values of U appear to be practically not renormalized in the range between n v H 2 and n = 1, otherwise the 
correlation effects destroy the magnetic ordering. 


C. Face-centered cubic lattice 

1. t' = 0 

The results for non-bipartite lattices should be substantially different, which is implied by Nagaoka theorem for the 
case of U —> oo££. In this section we treat the fee lattice. Since this lattice is not bipartite, the symmetry with respect 
to particle-hole transformation breaks down. Thus, for t' = 0 we show separately the HFA (Fig. flOh ) and SBA (Fig. 
[TUb ) results for the density range 0 < n < 2. The bare density of states is characterized by a divergent logarithmic 
singularity at the band top for t' = 0 (originating from the van Hove singularity line k = (0, Q , 27r), 0 < Q < 7 t/2 , V 
direction). The important difference between the fee lattice and those considered in the previous Sections is a good 
separation of the peculiar regions: n v h = 2 (van Hove singularity at the band top) and n = 1 (a very bad nesting) 
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Hartree-Fock Slave boson 



FIG. 10: (Color online) The same as in Fig. [3] for the face-centered cubic lattice Hubbard model with t' = 0 (a) HFA results; 
(b) SBA results. 


already at t' = 0. 

In the hole-doped half of the phase diagram we observe a strong suppression of any magnetic ordering by correlations 
up to n v m ~ 0.767 which corresponds to the L point of the Brillouin zone. At larger density we obtain a quantitative 
renormalization of the lower phase boundaries with the PM phase. At half-filling we observe a very narrow region of 
antiferromagnetism corresponding to the wave vector Q = (0, n, 2n) (a non-collinear AFM structure) which, however, 
appears to be unstable with respect to competition with another type of commensurate antiferromagnetism (wave 
vector Q = (0, 0, 2n)) at small hole doping and a vertically modulated FM phase (wave vector Q = (0,0, Q)) at 
small electron doping. Another correlation consequence is a moderate suppression of the FM ordering in favor of the 
diagonal phase (wave vector Q = (Q, Q , Q)). The FM ordering is saturated and appears to be stable well away from 
n = 2 only at U > 16 1. In close vicinity of the band top the case of a small number of carriers is realized. Therefore, 
the multiple scattering effects can substantially increase the critical value of U for the PM-FM transition: probably 
the most suitable way for the treatment of this problem is the T-matrix approach by Kanamori 60 . We believe that 
away from the band top (bottom) the slave-boson approximation yields a reasonable picture of magnetic ordering. 
The case of large U will be treated in detail elsewhere. 


Hartree-Fock Slave boson 



FIG. 11: (Color online) The same as in Fig. [3] for the face-centered cubic lattice Hubbard model with if = 0.3 1. 

To trace the impact of finite t' (which implies a finite DOS for the fee lattice due to the cutting of the logarithmic 
singularity), we calculated the ground state magnetic phase diagram in both HFA and SBA for if = 0.3f (Fig llll) and 
if = —0.31 (Fig. H2]). 
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Hartree-Fock Slave boson 



FIG. 12: (Color online) The same as in Fig. [3] for the face-centered cubic lattice Hubbard model with // = —0.3 1. 


2. t' = 0.3 1 

It is remarkable that for the case t' = 0.3 1 the correlation effects dramatically rebuild the phase diagram fully 
destroying any magnetic ordering in the vicinity of the band top and replacing the saturated FM ordering by the 
diagonal phase regions (Q = (Q, Q,Q) and (7r,7r,7r)) well away from it. The spiral and AFM phase regions are 
substantially narrowed far away from the value n = 1. There are 2 van Hove singularities in this case: n V Hi ~ 0.563 
(L point), n V H 2 ~ 1-87 (W point). Both singularities correspond to high symmetry points of the Brillouin zone, which 
seems to explain a rather small DOS for all the range of electron concentrations. Together with the absence of Fermi 
surface nesting this leads to a particularly strong suppression of magnetism by electron correlations as compared to 
all the diagrams studied. 


3. t' = -0.3 1 

On the other hand, at t' = — 0.3f there appears an additional saddle point with reduced symmetry n V H 3 ~ 1-9 
((0, a, a), cosa/2 = t/(2t' — t )) besides two other singularities u v hi ~ 1-12 (L point), n V H 2 ~ 1-33 (X point). The 
’vH3’ singularity leads to a relatively high DOS and consequently to low values of U c . As a result, renormalizations 
of the phase diagram in the region between n V H 3 and n = 1 become rather small. Above ro V H 3 the FM ordering is 
strongly suppressed, which dramatically disagrees with the Stoner criterion and HFA calculations, but agrees with 
the T-matrix approximation results^. However, when the number of carriers becomes equal to or less than n V H 3 we 
find a large saturated FM ordering region even at small U , which is suppressed by the spiral with Q = (0, 0, Q) in 
the vicinity of half-filling (this picture holds also within HFA). In general we conclude that the magnetic ordering 
in the electron-doped half of the phase diagram is only quantitatively affected by the correlations: even the Stoner 
criterion works fine in this region. In contrast to the other lattices the van Hove scenario does not determine wave 
vector in vicinity of n V H 3 ’s in general case, the high DOS factor prevails in the electron-doped half and causes a 
tendency to structures with small wave vectors. In the hole-doped half of the phase diagram we find the typical for 
SB A suppression of any magnetic order far away from n = 1. What is unusual in the hole-doped half, is the presence 
of wide regions of PS between the PM and different spiral magnetic states. 


D. Analysis of bosonic fields: comparison of Hartree—Fock and slave-boson approximations 

To illustrate the correlation effects impact in detail, we show the concentration dependence of site occupancies for 
the forced PM state and the magnetically ordered state within both HFA (given by Eq. (l25jl ) and SBA (Eq. (l22l) ). 
Figure [13] depicts the behavior of bosonic variables depending on the electron concentration n for the sc lattice with 
t' = 0.3f, U = 9t for the PM phase (a) and the spiral (antiferromagnetic) ground state (b). We present d 2 for n < 1 and 
e 2 for n > 1, since doubles and holes play the role of current carriers. Below we discuss for brevity the case of n < 1, 
analogous conclusions hold for n > 1 with the replacement d —> e. In the PM case d 2 highly increases as n approaches 
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FIG. 13: (Color online) The density dependence of the physical quantities for simple cubic lattice with t' = 0.3 1 and U = 9t for 
(a) paramagnetic state and (b) magnetically ordered ground state. Upper panel: The magnetic moment m calculated within 
SBA (HFA) is shown by black solid (dashed) line on left axis. Filled (empty) circle shows the magnetization for AFM phase 
at n = 1 within SBA (HFA). The energy difference between HFA and SBA A E = Ebf — Esb is shown by purple dash-dotted 
line on right axis. The cross shows this energy difference for AFM phase at n = 1. Middle panel: The renormalization factor 
z 2 (z|) is shown by blue solid (dashed) line; the circle shows it for AFM phase at n = 1. Qx,z/n is shown by green solid line. 
Lower panel: Occupancies of the double carrier (d 2 for n < 1 and e 2 for n > 1) calculated within SBA (HFA) are shown by 
solid (dashed) red line; filled (empty) circle shows it for AFM phase at n = 1 calculated within SBA (HFA). 


the half-filling within HFA, which strongly disagrees with its suppression (and the corresponding bandwidth narrowing 
z 2 ) within SBA. The latter results in the metal insulator transition for n — 1 at a sufficiently large U/l£2- when the 
spectral weight z 2 vanishes. This treatment of the metal-insulator transition is physically close to the Brinkman Rice 
one^i. Therefore, within SBA, in the PM state electrons avoid double site occupancy, which allows their Coulomb 
energy to be lowered (the price is an increase in kinetic energy due to the band narrowing): in bosonic language this 
means a much more favorable way of redistribution between p 2 and d 2 than in HFA (cZsba ^hfa)- 

A very different picture can be found for the magnetically ordered case, where the system has an additional degree 
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of freedom (redistribution between pi and pi). We found that already within HFA the Coulomb energy turns out to 
be substantially reduced due to the lowering of d 2 , and the distinction from the SB A in this sense is small. In this 
case one finds that the main correlation effect is the redistribution between p 2 and p 2 which results in: (i) a partial 
magnetization suppression, (ii) a moderate bandwidth reduction z 2 due to the Pauli principle. The general conclusion 
is that the correlation impact on the magnetically ordered state appears to be much more weaker than on the PM 
phase. Therefore, the main “correlation” effect (in a broad sense) is the magnetism formation itself. The data on 
magnetic moment and z-factors of Fig. 113b demonstrate the jump of magnetic moment near first order transition 
(Am ss 0.2 for n < 1 and Am ~ 0.4 for n > 1). The renormalizations of effective masses given by z~ 2 are considerably 
spin-dependent. In particular, at large U in the ferromagnetic region we have z 2 ~ 1 and z 2 ~ 1 — n = e 2 . In fact, in 
this limit the spin-down states have non-quasiparticle nature^. 

For the half-filling e 2 and d 2 are equal, being small (Fig. I13b h which implies that the spectra renormalization 
factor z 2 tends to 1 at n = 1. The increase of U causes the change of the gap nature from the Slater to Hubbard one, 
the latter being retained in the paramagnetic state. In particular, the difference between mass renormalization within 
ordered and paramagnetic case was discussed in Ref. [66(. However, near half-filling HF works well for the electron and 
spin-wave spectrum at U |t|, yielding exactly atomic limit 6 '. While in the ordered state the mass enhancement, 
1 /z 2 , is smallest at n = 1, at doping this quantity rapidly grows as one can see from Fig. H3b. This explains why the 
SBA results are very similar to the HFA ones in the vicinity of half-filling: the sequences of the magnetic states on 
the diagrams do not change with increasing U, although their regions shift to larger C/’s in SBA (Fig. [T5bh 

Upper panels of Fig. Q3] demonstrates also the energy difference A E = Ehf — E$b between HFA and SBA. A E 
for PM phase reaches 0.63C and is much larger than A E for magnetic state, which reaches 0.23/ only. This illustrates 
also the conclusion that correlation effects lower the energy of PM phase much stronger than the energy of magnetic 
states. The energy lowering at half-filling is maximal for PM phase and minimal for spiral state, which agrees with 
the dependencies of z-factors. 

We note that for other lattices the results for bosonic fields and the energy difference behavior are qualitatively the 
same. 


VI. CONCLUSIONS 

We have investigated the ground state magnetic phase diagram including spiral ordering formation and phase 
separation within both Hartree-Fock and slave-boson approaches. The slave-boson approximation enables one to take 
account of the correlations in terms of a few parameters renormalizing the electron spectrum. The main advantage of 
SBA is a correct qualitative description of the paramagnetic phase: the energy becomes lower as compared with HFA, 
so that the gain in energy of the magnetic phases is now substantially smaller. As a result, the paramagnetic region 
occurs in the ground state within the SB method at large U. Being technically simpler in comparison with DMFT, SB 
mean-field approximation gives a transparent picture of electron spectrum, but misses dynamical self-energy effects. 
However these effects can be in principle taken into account within SBA beyond mean-field level. 

Spiral configurations within DMFT for the square lattice with t' = 0 were treated in Refs. [36ll37i . The sequence of 
phase transitions for U = 8/ is AF—>• AF+(Q, Q) -» (Q, Q) —> (Q, Q) + ( Q , n) —> (Q, ir) —»• PM. This is in agreement 
with out results, although the corresponding U value is somewhat larger. According to Ref. [3(1 magnetic ordering 
vanishes for n < n c = 0.35, whereas our critical value n c is considerably larger, about 0.6. Recent inhomogeneous 
DMFT calculations^ enabled one to describe formation of spin density wave and stripe states in the 2D case. The 
results are in a qualitative agreement with our data (see Fig. [2J. A quantitative difference of the critical concentration 
(n c = 0.8 instead of our value 0.6 for U > 10/) may be due to that we consider spirals rather than spin density waves. 
In Ref. [H a comparison of spiral and stripe states was performed within the SB method for different n and t'/t. The 
spirals turn out to be more energetically favorable for large t'/t and n = 7/8. Besides that, stripe states may are 
spatially homogeneous and loose in long-range Coulomb interaction energy in comparison with spiral states. 

Our results and their comparison with DMFT enable one to estimate the role of dynamical corrections to the 
electronic spectrum, since the SBA treats the nature of correlated electron states only on the static level. 

The results obtained for both 2D and 3D cases can be applied to an analysis of the properties and phenomena 
in d-metals and their compounds, especially with fee and bcc structures. The correlation effects considered can be 
important, e.g. for the metal-insulator (Mott) transition (this problem was treated in Ref. [34] in the Hartree-Fock 
approximation including the spiral states). 

We see that the correlation effects lead to a strong suppression of the regions of the existence of the magnetic phases. 
At the same time, there remain the first-order transitions and noticeable regions of the phase separation between the 
magnetically ordered states. The correlation effects near half-filling change only slightly the Hartree-Fock results, so 
that at small Z'/Z the sequences of magnetic states do not change with increasing U. 
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Both SBA and HFA yield qualitatively (but not quantitatively) similar results near half-filling, including phase 
separation. The main advantage of SB method with respect to HFA is correct account of the doubly occupied sites 
onto the one-electron spectrum which manifests itself through the 2 -factors. Owing to them PS regions become more 
narrow within SB method. 

Our results confirm the importance of the van Hove and nesting scenarios in stabilization of spiral magnetic 
structures. The divergent van Hove singularities which are characteristic for 2D lattices, as well as for bcc lattice with 
t' = 0, prevent the suppression of magnetism by the electron correlations. The role of spectrum renormalization within 
SBA in comparison to HFA is weaker when the physics is governed by van Hove singularities. Far from half-filling, 
the van Hove singularities stabilize the ferromagnetic order, especially in the 2D case. The correlation effects are 
shown to be negligible for the FM state in the case of a large density of states, and for the spiral states provided that 
the Fermi surface nesting with corresponding Q is present. The most important among the bounded singularities are 
those produced by the saddle points which do not correspond to high symmetry points of the Brillouin zone. They 
favor stabilization of incommensurate spiral states and usually lead to an increase in the density of states. The value 
of wave vector Q is equal to a vector connecting two of the saddle points in these cases. When the system is far from 
half-filling and the Fermi level is far from van Hove singularities, the magnetic state cannot form at any U. 

A recent series of paper s 68 ’ 69 concerns the formation of spiral magnetic order in a specific many-band model (two- 
band Rice model 7 -) with imperfect nesting between the hole and electron Fermi surface pockets. The general result of 
these investigations is that there exists a first order transition between commensurate and incommensurate magnetic 
states which invokes a phase separation between them. Thus the phase separation is a more general phenomenon and 
the existence of two nested Fermi surfaces is not necessary. 
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